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Abstract. The annihilation of neutralino dark matter in the Galactic Center (GC) 
may result in radio signals that can be used to detect or constrain the dark matter 
halo density profile or dark matter particle properties. At the Galactic Center, the 
accretion flow onto the central Black Hole (BH) sustains strong magnetic fields that 
can induce synchrotron emission by electrons and positrons generated in neutralino 
annihilations during advection onto the BH. Here we reanalyze the radiative processes 
relevant for the neutralino annihilation signal at the GC, with realistic assumptions 
about the accretion flow and its magnetic properties. We find that neglecting these 
effects, as done in previous papers, leads to the incorrent electron and photon spectra. 
We find that the magnetic fields associated with the flow are significantly stronger than 
previously estimated. We derive the appropriate equilibrium distribution of electrons 
and positron and the resulting radiation, considering adiabatic compression in the 
accretion flow, inverse Compton scattering off synchrotron photons (synchrotron sclf- 
Compton scattering), and synchrotron self- absorption of the emitted radiation. We 
derive the signal for a Navarro-Frenk- White (NFW) dark matter halo profile and 
a NFW profile with a dark matter spike due to the central BH. We find that the 
observed radio emission from the GC is inconsistent with the scenario in which a spiky 
distribution of neutralinos is present. We discuss several important differences between 
our calculations and those previously presented in the literature. 



1. Introduction 

The combination of cosmic microwave background anisotropy studies with observations 
of distant type la supernovae and measurements of the large scale structure of the 
universe reveal that the energy density in our universe is dominated by dark energy 
(~ 70%) followed by a significant contribution from matter (~ 27%). In addition, these 
studies show that the baryonic content of the universe is limited to less than ~ 4% in 
agreement with big bang nucleosynthesis predictions. More specifically, recent WMAP 
results give the mean baryonic density in the universe to be within 0.044 ±0.004 of the 
critical density while the matter density is 0.27 ±0.04 of critical The remaining ~ 
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0.23 of the matter density, i.e., the bulk of the matter in the universe, is beheved to be 
a yet undiscovered form of dark matter. 

Weakly interacting massive particles (WIMPs) are natural candidates for the 
dark matter j2]. Particles with masses around ~ 100 GeV that interact only weakly 
have freeze-out densities in the required range. In addition, particle physics models 
that invoke supersymmetry generate a number of plausible WIMP candidates. In 
supersymmetric extensions of the standard model, the lightest supersymmetric particle 
may be stable due to conservation of R-parity enabling their survival to the present. 
In addition to massive and weakly interacting, dark matter particles are expected to 
be neutral. A class of neutral lightest supersymmetric particles is a combination of 
gauginos and higgsinos, named the neutralino often represented by x- 

A number of experiments are presently searching for the neutralino (see, e.g., jSj). 
Accelerator experiments have placed lower limits on the neutralino mass of mj^^ZI 
GeV 13] . The neutralino mass in the minimal supersymmetric extension of the standard 
model that can explain the dark matter density is also bounded from above, 
TeV. In some constrained supersymmetric models m^^SOO GeV due to the cosmo logical 
constraints from WMAP [Hj. Direct searches for neutralinos that cross the Earth are 
presently reaching the required sensitivity to probe the relevant range in parameter 
space of supersymmetric models. 

A good complement to accelerator and direct detection methods are indirect 
searches based on the emission from neutralino self-annihilation in astrophysical 
environments |EllZllHlinilIIIlinilI2llinilI31linilin!. Neutralinos self-annihilate at a 
rate proportional to the square of the neutralino density. Thus, the highest density 
dark matter regions are the best candidates for indirect searches. In fact, the GC region 
may potentially be so dense that all neutralino models would be ruled out 0. This 
strong constraint arises in models where the super-massive black hole (BH) at the GC 
induces a strong dark matter density peak called the spike. The existence of such a 
spike is strongly dependent on the formation history of the Galactic Center BH If 
the BH formed adiabatically a spike would be present, while a history of major mergers 
would most likely not allow the survival of a spike. 

In contrast to the uncertain presence of a central spike in the dark matter 
distribution, the central BH is known to induce an accretion flow of baryonic matter 
around its event horizon. The accretion flow carries magnetic fields, possibly amplified 
to near equipartition values due to the strong compression. The distribution of electrons 
and positrons (hereafter both called electrons) produced by neutralino annihilation at 
the GC would also be compressed toward the BH radiating through synchrotron and 
inverse Compton scattering off the local photon background. 

By considering the injection of electrons, combined with radiative losses and 
adiabatic compression, we find the equilibrium spatial and spectral electron distribution 
and derive the expected radiation signal. We find that the synchrotron emission 
of electrons from neutralino annihilation ranges from radio and microwave energies 
up to the optical, in the central more magnetized region of the accretion flow. At 
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low frequencies, synchrotron self-absorption (SSA) reduces the amount of radiation 
transmitted outwards. The resulting signal is stronger than the observed emission 
in the 10 to 10^ GHz range for the case of a spiky dark matter profile while for a 
pure NFW profile the emission is below the observed level. Here, we do not attempt 
to explore all possible spike profiles in the GC region, that might originate from the 
adiabatic compression of different types of initial dark matter profiles. Instead, we limit 
ourselves to the case of a spike arising from adiabatic compression of a NFW profile 
in the gravitational field of the massive BH at the GC (hereafter, the spike case.) in 
contrast to the pure NFW profile. 

Previous calculations of the type presented here were either incomplete or had 
unphysical values for some of the parameters. In particular, calculations presented in 
[Oj and [l2j assume a magnetic field that does not correspond to the equipartition field 
for the accretion flow around the black hole. We show here that if this incorrect field is 
adopted, then the main processes responsible for reaching the equilibrium of relativistic 
electrons are synchrotron self-compton scattering (SSC) and adiabatic compression of 
the particles advected with the accretion flow, processes which were both ignored in 
previous calculations. On the other hand, if the correct magnetic field is adopted, 
the equilibrium distribution of the radiating electrons is determined by advection and 
synchrotron emission. 

The paper is organized as follows: in ^ we describe the spatial distribution of 
dark matter in the GC region. In ^ we describe the accretion flow around the 
BH and its magnetic field. The injection of electrons and positrons from neutralino 
annihilation is considered in 2] In ^ we solve the transport equation for electron- 
positron pairs from neutralino annihilation in the presence of energy losses and adiabatic 
compression. Synchrotron self-absorption is discussed in ^ Results of our calculations 
and comparison with observations are shown in ^ and comparison with previous work 
is made in §H1 We conclude in ^ 



2. The spatial distribution of dark matter in the Galactic Center 

The spatial distribution of dark matter in halos is still a matter of much debate (see, 
e.g., [IS])- Numerical simulations suggest that coUisionless dark matter forms cuspy 
halos while some observations argue for a fiat inner density profile [20] ■ The standard 
numerical dark matter halo is the NFW profile [21] which is expected to be universal. 
However, recent simulations have found more cuspy halos [22] as well as shallower profiles 
[2H]- At present, it is not clear that there is a universal dark matter halo profile, but the 
NFW case seems to represent well the range of plausible profiles. Therefore, we assume 
that the NFW profile describes well the dark matter in our Galaxy. 
The NFW profile can be parametrized as 

pW = /° ,2 - (1) 
— 1 + — 

where Rc is the core radius and po is a normahzation constant. The two free parameters 
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can be determined from the total mass of the Galaxy (that we take to be 2 x lO^^M©) 
and the velocity dispersion at the solar location which gives = Pdm{Rq) = 
6.5 X 10-25g/cm^ for Rq = 8.5 kpc. 

There is now growing evidence for the presence of a supermassive black hole at 
the Galactic Center, with mass ~ 2 x lO^M©. In fact most galaxies seem to harbor 
central black holes with comparable or even larger masses. The presence of a BH can 
steepen the density profile of dark matter by transforming the cusp at the GC into a 
spike of dark matter ^24]. The density profile of the spike region, where the gravitational 
potential is dominated by the BH is described by 



p'w = «r-M— ^ P0^?w — • (2) 



Here, 6 is the slope of the density profile of dark matter in the inner region (5 = 1 for a 
NFW profile), and 6sp = (9 — 26) /(4 — 5). The coefficients as and g{r) can be calculated 
numerically as explained in detail in [21]. It is possible to identify a spike radius Rsp 
where the spike density profile given in Eq. matches the NFW dark matter profile. 
In other words, at Rgp the gravitational potential is no longer dominated by the central 
BH. 

Neutralino annihilations affect the density profile in the spike by generating a 
fiattening where the annihilation time becomes smaller than the age of the BH. This 
fiattening occurs at the position 

Pcore 7 ^ 7 ) K,"^) 
\<yV)anntBH 

where {crv)ann is the thermally averaged cross section times velocity for neutralino 
annihilation. Cosmological arguments give estimates of (crf)ann~3 x lO^^^cm^s"^. The 
effect of annihilations on the spike density profile can be written as 

^) = / / N , ' (4) 

Psp{r) + pcore 

which accounts for the fiattening in the central region. 

Several dynamical effects may weaken or destroy the spike in the GC j21 El ITH] , 
depending on the history of formation of the central BH. If the spike is not formed 
or gets destroyed, the central region should be described by the cuspy profile as in the 
NFW case. Here we consider both cases and show that the observed emission is stronger 
than the predictions for a NFW cusp, while the spike generates a signal well above the 
observed data. 



3. The accretion flow and its magnetic field 

We model the accretion fiow of gas onto the black hole at the center of our Galaxy 
following a simple approach described in |25j. More detailed models of the accretion 
flow around the BH lead to corrections which are negligible when compared to the 
uncertainties in the dark matter distribution. In the model we adopt, the BH accretes 
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its fuel from a nearby molecular cloud, located at about 0.01 pc from the BH. The 
accretion is assumed to be spherically symmetric Bondi accretion with a rate of mass 
accretion of M = 10^^M22 g s~^. The accretion onto the BH occurs with a velocity 
around the free-fall velocity, such that 



v{r) = PGMbh/t = c ( ^ ) (5) 



where Rg = 2GMbh/c^ = 7.4 x 10^^{Mbh/'2-5 x IO^Mq) cm is the gravitational radius 
of the BH and Mbh is the BH mass. Therefore, 

t;(r) = 1.46 X lOM \ cm s-\ (6) 

^ ' \2.5 X lOeMo/ l^O.Olpc/ ^ ^ 

Mass conservation then gives the following density profile: 

M _ M f r 
ATTr^v{r) AttR'^c \Rg 

such that 



pir) = 5.6 X 10-'^M22 g cm-3 . (8) 

Following [23 5 we assume that the magnetic field in the accretion flow achieves 
its equipartition value with the kinetic pressure, namely pv'^/2 = i?(r)^/8'7r. With this 
assumption. 

It is believed that magnetic fields in the accretion flow will in general reach the 
equipartition values described in Eq. 0. However, smaller fields may be reached if the 
equipartition is prevented somehow. In fact, previous authors P [121 have assumed a 
lower field given by 

/ \-5/4 / \-5/4 

Bio^{r) = IpG [ ^ j = 3 X 10^ pG [j^] - 10"' , (10) 

although this field was actually referred to as the equipartition field. In what follows, 
we assume the equipartition field as calculated in Eq. and derive the signal from 
neutralino annihilation. In ^ we address the case of assuming a lower filed such as in 
Eq. (Hm). 



4. Neutralino annihilations: the injection of electrons and positrons 

Neutralino annihilations can result in many different final states (see, e.g. ^Hj and 
references therein). Given the large uncertainties in the dark matter distribution, we 
chose to make the simplifying assumption that the annihilation channels are dominated 
by quarks and antiquarks, which in turn generate mainly pions: 

X + X-" (1 + ^-^ pions. (11) 
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Approximately 1/3 of the pions are neutral and promptly decay to gamma rays. About 
2/3 of the pions are charged and result in the production of electrons, positrons and 
neutrinos. Although much attention has been devoted in the past literature to the 
gamma rays from vr*^ decays, electrons and positrons can result in copious production of 
nonthermal radiation in the presence of magnetic fields and a photon background. 

We simplify the annihilation treatment further by assuming that the spectrum of 
charged pions in a single neutralino annihilation is described by the Hill |26j quark 
fragmentation: 

W^{E^) = -^—x-^'\l~x)\ (12) 

where x = ET^/m^. In Eq. fll2|) . a 2/3 factor accounts for the fraction of charged 
pions compared to neutral (1/3) in each neutralino annihilation event, and a factor 2 
accounts for two jets, each with energy m^, the neutralino mass. With these simplifying 
assumptions, the spectrum of generated electrons appears to possess the main features 
found with a full treatment of the neutralino annihilations (see |13j for a description 
of the full calculation), although the latter should be used for a more quantitative 
description of the particle physics aspects of the problem. 

The spectrum of electrons (and positrons) from the vr^ decays is calculated by 
convoluting the spectrum of pions and muons. For relativistic electrons the electron 
spectrum reads 



We{E, 



dE„ 



max{i?e,mmu) 
2 

mi 



dE^W^{E^) ■ 
1 dueiEe, E^, E^) 



- ^|El - ml 



dE, 



(13) 



where, neglecting the muon polarization, we get 



dUejEe, E^, E^) 
dK 



E^P 



5 _ 

6 

(1-/3)2 



3^2 _|_ 2 3 
2^ 3 



if i-f<.<l 

if 0<6<i^, 



(14) 



with e = Y+/3^- Here P is the pion speed and E^, and E^^ are the total energies of 
electrons and muons respectively. The two limits of integration i?™™(£'^) and E^'^^lE^) 
can be derived by inverting the following equations: 
E^ 



2mi 



E^ 



m. 



2m£ 



iii + P) + Kil - P) 



mUl-P)+ml{l + P) 



(15) 



(16) 



Finally, the injection of new electrons at the distance r from the BH at energy E 
can be written as follows: 

2 



Q{E,r) 



1 / PDMjr) 

2 V mv 



W,{E){av), 



(17) 
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where the density of dark matter has the profile pDhiij") discussed in ^ with or without 
a spike. In this equation, the factor 1/2 is introduced to avoid double counting of the 
annihilating neutralinos, as pointed out in 



5. The equilibrium distribution of electrons near the BH 

The spectrum of particles at a position r in the accretion flow around the BH is the 
result of the injection of newly produced electrons at the same position, radiative losses 
of these electrons, and the adiabatic compression that may enhance their momentum as 
they move inward. Here we neglect spatial diffusion, which occurs on larger time scales. 
The transport equation including all these effects can be written as follows: 

- [r''v{r)\ + 1^ [p'p{r,p)f] = Q{r,p), (18) 

where f{r,p) is the equilibrium distribution function of electrons injected according to 
Q{r,p) that lose energy radiatively as described by the function p{r,p) = dp{r,p)/dt. 
Here v{r) = —c{r/Rg)^^^'^ is the inflow velocity. The equation can be solved analytically 
if the electrons remain relativistic everywhere in the fluid, which is a good approximation 
for the range of frequencies of the radiation that we are interested in. 
In Eq. ()18|). the term 

1 1 (9 r 1 

Pad = --pVw(r) = -^P^ v{r) (19) 
6 or l J 

describes the rate of change of momentum of a particle at the position r due to adiabatic 
compression in the accretion flow. The rate of adiabatic momentum enhancement should 
be compared with the rate of losses due to synchrotron emission: 

4 B'^(r) 

Psyn{r,p) = -(JT^^l^ (20) 

where ax is the Thomson cross section and 7 is the Lorentz factor of the electron. The 
magnetic field Beg{r) depends on r as described in Eq. 0. Comparing the absolute 
values of the adiabatic compression rate and the rate of synchrotron losses, we find that 
the effect of adiabatic compression dominates at 

r 

- > 32.7 7 . 

Rg 

For this estimate we used the fiducial value of the BH mass (2.5 x IO^Mq) and the 
accretion rate (M = lO^^g s^^). A graphical comparison of the time scales for adiabatic 
compression (solid lines) and synchrotron losses (dashed lines) is plotted in Fig. Qfor 
7 = 10 and 7 = 1000. 

It is clear from Fig. ^that for Lorentz factors 7 < 1000, there is a region (large radii) 
for which the particles injected at r are advected inward faster than they lose energy 
through synchrotron emission. If instead of using Beq{r) from Eq. we adopted 
Biow{r) from Eq. (fTIH) as in [HI US], the relation for the dominance of the adiabatic 
compression would become 

^ > 2.3 X 10"^ 7 . 

Rg 
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In this case, the adiabatic compression (ignored in previous calculations) becomes more 
important than synchrotron losses at most energies. The final spectrum would then 
differ substantially from that obtained in [^1 , as we discuss in ^ 

In order to solve the transport equation, Eq. (|T8|l . we discuss different loss processes 
using the equipartition field Beq. We first consider synchrotron losses (psyn), followed 
by inverse Compton losses ijpics) and synchrotron self-compton scattering ijpssc)- 

The rate of synchrotron losses from Eq. (pn|) is given by: 

Psyn{r,p) c = 1.6 X 10 L Q^^^ 1 7^ erg s ^ . 

In the Thomson regime, losses due to Inverse Compton Scattering (ICS) off a photon 
background with energy density Uph has the following form 

4 

Vicsir.p) = -aTUphir)-i'^ ■ (21) 

ICS dominates synchrotron losses only if > SirUph- 

Assuming that Uph is independent of r, i.e., a fixed photon background, ICS becomes 
important at large radii (~ O.Olpc) and only if f/p/^^lO^eV cm~^. Such strong photon 
background is unlikely to be present at the GC region. For comparison, the CMB 
radiation has Ucmb ~ 0.25eV cm~'^ while the optical background has Uopt ~ leV cm~^. 
If ICS off a fixed background is not dominant at large radii, it becomes even less 
important at small radii when compared to synchrotron losses. Consequently, we safely 
neglect the role of ICS off photons of fixed photon backgrounds. 
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The electrons, radiating in the strong magnetic field near the BH generate a photon 
background that can become quite intense. The rate of losses due to ICS of electrons 
off the photons generated through synchrotron emission by the same electrons is 



Pssc{r,p) = gC^rt^pT('^)7^ 



(22) 



where the photon energy density generally depends on the radius r. The photon density 
Up[^{r) is a nonlinear function of the distribution f{r,p). In other words, the term p in 
Eq. ()18|) depends on f{r,p) when synchrotron self-Compton scattering is included. If 
this contribution cannot be neglected, an analytical solution of the transport equation 
becomes unattainable. 

Given the distribution function f{r,p), one can calculate the synchrotron emissivity 
j(i^, r) (energy per unit volume, per unit frequency, per unit time). The photon energy 
density at the position r is then proportional to the integration over all lines of sight of 
the emissivity, with the possible synchrotron self-absorption taken into account at each 
frequency. However, the distribution function f{r,p) is not known a priori, and the 
problem becomes intrinsically nonlinear. The approach that we follow here is to start 
with neglecting synchrotron self-Compton scattering and check a posteriori whether the 
assumption is correct in the situation at hand. 

An analytic solution of the transport equation can be derived when p{r,p) is 
dominated by synchrotron losses as in Eq. (^0]) . In this case, the equation admits 
the following analytical solution: 

-2 



f{r,p) 




Rac 



dRi 



R. 



5/2 



Rn 



'pinj[p, r, R 



P 



Qi^Rinj ; Pinj \P: Rinj^ ) • 



(23) 



We rewrite the synchrotron losses in the form: 



-5/2 



Psyn = ^0 ( 1 



where = aTB^/QT^ijncY . The function Pinj\p,'r, Rinj] corresponds to the injection 
momentum of an electron injected at the position Rinj that arrives at the position r 
with momentum p. This injection momentum can be obtained by inverting, with respect 
to Pinj the solution of the equation of motion of the electron, in the presence of adiabatic 
compression and radiative losses: 

Pi^) ■ (24) 



dp 
dr 




p'- 



2Ra 



^9 J ^-"-9 \^^9 / 

The solution of this equation, with initial condition p[r 



p[r,pinj,R, 



Pinj 



2koRl 
3c 



Pinj 



Rinj J Pinj: Rinj^ Pinj is 
\ 3/21 / \ 1/2' 



1 



R 



in] , 



Ri 



(25) 
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Figure 2. Electron density per unit energy as a function of energy (solid lines) at 
r = 10^ Rg (upper curve) and r = 10^ Rg (lower curve). Superimposed (dashed curves), 
wc plot the function n{E, r) obtained without the effect of advcction. The left panel 
illustrates the spike case, while the right panel applies to the NFW case. 



In the absence of synchrotron energy losses (when kg — 0), particle momenta only change 
due to adiabatic compression, and the momentum of a particle changes according with 
the well known p — Pinj{r/Rinj)~^^'^, vahd for the case of free fall. 

The joint effect of the energy gain due to the adiabatic compression and the 
energy losses due to synchrotron emission generates a new energy scale in the system. 
Introducing 

3rc 97rr(mc)^ 
we can rewrite p[r,Pinj, Rinj] as: 

P\j'iPinji Rinj\ Pinj 

Prom this expression it is clear that, at any fixed position r, adiabatic compression 
dominates over synchrotron losses if the injection momentum is appreciably lower than 
Pm- In this case, the electron energy increases while the electron moves inward, until the 
rate of synchrotron losses becomes important. The opposite happens when the electrons 
are injected at momenta larger than p^, since synchrotron losses are important from 
the time of injection. 

The momentum pm can be interpreted as the momentum where the two competing 
effects of adiabatic heating and synchrotron losses balance each other. Thus, particles 
accumulate at momentum pm- This phenomenon depends on the distance from the GC: 
at large distances from the BH the momentum p^, which scales linearly with radius, is 
large and the local rate of injected electrons is low, therefore, the accumulation is small. 
At small distances the accumulations at Pm becomes more evident. 

We define now the equilibrium electron spectrum n{E, r) at the position r, which 




(27) 
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is related to f{E,r) through the relation 

n{E, r)dE = 47rp^/(p, r)dp. 

In Fig. 121 we plot n{E,r) as a function of the electron energy at two different radii, 
r = 10^ Rg (upper solid curve) and r = lO'^Rg (lower solid curve). The dashed lines 
illustrate the solution of the transport equation when adiabatic compression is switched 
off and only synchrotron losses are included (as in ^2])- The dark matter density profile 
is assumed to be characterized by a spike at the center (2^ in the left panel of Fig. 
121 and by an NFW density profile |21j in the right panel of Fig. |21 The neutralino 
mass assumed in Fig. |21was 100 GeV, while the annihilation cross-section was chosen 
to be 10~^^cm^s~^. The accumulation effect described above manifests itself through 
the appearance of the spiky feature at the momentum pm- The accumulation is less 
pronounced at large radii, as expected. 

We conclude this section by addressing the issue of the synchrotron self- Compton 
scattering. As explained above, this effect cannot be accounted for in an analytical 
approach to the transport equation, since it is intrinsically nonlinear. Instead, we check 
a posteriori if neglecting SSC was a good assumption. Given the symmetry of the BH 
region, the photon energy density as a function of r can be written as: 

1 



Tjsyn, 



(28) 



where j(z/, r) is the synchrotron emissivity. 

This energy density can now be compared with the magnetic energy density at the 
same location, i?^(r)/87r. The results are plotted in Fig. IHl the solid line represents the 
photon energy density Upl"'{r), while the magnetic energy density is plotted as a dashed 
line. The calculations are carried out for a dark matter density profile with the spike at 
the center and a neutralino mass of 100 GeV. The magnetic energy always dominates 
over the photon energy, although at large distances from the BH the separation between 
the two curves reduces to about one order of magnitude. The curves in Fig. 01 are 
obtained without taking into account the synchrotron self-absorption (see next section) 
at low frequencies, therefore the calculated photon energy density should be considered 
as an upper limit. On the basis of this argument, we can conclude that neglecting SSC 
was a reasonable approximation for the scenario at hand. 



6. Synchrotron Self- Absorption 

Synchrotron radiation can be reabsorbed by the radiating electrons when the system 
is sufficiently compact. This phenomenon of synchrotron self-absorption (SSA) is 
particularly effective at low frequencies. The radiation transfer across the GC region 
can be described by P7j : 

= -x{v, s)Iiu, s) + ^jiu, s) , (29) 
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Figure 3. Comparison between the energy density of the synchrotron emitted photons 
Uph (soHd Une) and the energy density associated to the magnetic field i?^/87r (dashed 
hne). 



where /(z/, s) is the radiation intensity at frequency v along a fixed hne of sight, s is the 
hnear coordinate along the line of sight, xi''^) s) is the absorption coefficient and s) 
is the emissivity. 

For sufficiently high frequencies, where the absorption coefficient is small, the 
radiation intensity is simply 

^ ^dsj{u,s). (30) 



/(z/,s) 



477 



The absorption coefficient can be written in the form: 



X{^, s) 



J dEP{i^,E)E^' 



_d_ 

dE 



n{E, r(s)) 
^2 



(31) 



where P(z/, E) is the synchrotron power radiated per unit frequency by one electron of 
energy E, and n{E, r(s)) is the electron density per unit energy discussed in the previous 
section, evaluated at the position r(s). Integrating equation (pTTjl by parts, using the 
vanishing electron density at E = we obtain: 

'■ ^n{E,r) d 



where 



P{u,E) 



dE- 



E^ dE 



E^P{p,E) 



mc^ 



ME). 



Here K^i^ is the modified Bessel function of order 5/3 and 



( E 



,mc^ 



(32) 



(33) 



(34) 
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Figure 4. Optical depth as a function of the photon frequency in two cases: 1) 
advection plus synchrotron losses (solid line); 2) synchrotron losses only (dashed line). 
The calculation was carried out for a line of sight tangent to the event horizon of the 
black hole. 



The absorption coefficient is therefore 

The effect of the SSA can be assessed by evaluating the opacity to the absorption 
process at fixed frequency defined as 

Tyih) = J c/sxhr(6, s)], (36) 

where b is the distance of the hne of sight from the hne crossing the GC. 

Clearly SSA becomes important when T^{b) approaches or exceeds unity. When 
Tuib) <ti 1, SSA can be neglected and the radiation intensity can be simply evaluated 
as in Eq. Whether the SSA is important or not depends also on the line of 

sight considered: in general the GC can become opaque to SSA at different frequencies 
for different lines of sight. The observed radiation is then the result of a line of sight 
integration of the radiation intensity, as described below. 

In Fig. m we plot the opacity Tj, for a line of sight tangent to the event horizon 
of the BH (solid line), namely with b = Rg. In this figure, we chose the dark matter 
density with a spike, = 100 GeV and {crv)ann = 10^^^ cm^/s. The dashed line 
in Fig. 0] represents the opacity to SSA if only synchrotron losses are included while 
adiabatic compression is neglected (as in P2])- For the correct calculation in which both 
synchrotron losses and advective compression are included, the SSA becomes important 
for frequencies lO^'' — 10^^ Hz. When only synchrotron losses are considered, SSA 
becomes important only at much lower frequencies. 



Neutralino Annihilation at the Galactic Center Revisited 



14 




Figure 5. Emitted luminosity in two cases: 1) advection and synchrotron losses (solid 
line); 2) only synchrotron losses (dashed line). The computation is performed with 
rrix — 100 GeV and {av)ann = 10^^^ cm^/s. The left panel represents the Spike case 
while right panel the NFW case. 



The dependence of the SSA on neutrahno parameters scales with n{E, r) oc 
{av)ann rn~^l'^ . At fixed {av)ann, the electron density drops when increases, and 
SSA decreases as well. Finally, in the case of an NFW density profile without a spike, 
the rate of electron injection remains such that the effect of SSA can always be neglected 
at the frequencies of interest. 

7. Results 

In this section, we present our results in terms of the synchrotron luminosity produced 
by relativistic electrons from neutralino annihilations in the GC. We considered both 
cases of a spiky density profile and a NFW density profile. In order to compute the 
radiation intensity 1^, we solved Eq. ^I^i in the case of the spiky density profile in 
which, as discussed in the previous section, the SSA effect is relevant. In the case of the 
NFW density profile we simply used equation pUjl . Integrating the radiation intensity 
over the sky (namely, over all lines of sight to the observer), one obtains the luminosity: 

L^ = 47r J db27rbdbh{b) . (37) 

In Fig. "we fix = 100 GeV and {crv)ann = 10"^'' cm^/s, and show the luminosity for 
the spiky case (left panel) and the NFW density profile (right panel). The continuous 
line shows the case in which advection is taken into account, while the dashed line shows 
the case in which only the synchrotron energy losses are accounted for. 

A comparison of the dashed and continuous curves in Fig. shows that the 
combined effect of advection and synchrotron losses appears at frequencies up to 10^^ Hz. 
In this frequency range, the electron density is much higher when advection is included 
as compared to the pure synchrotron case. This effect produces an enhancement of 
the emitted luminosity of about one order of magnitude. In addition, we can see that 
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Figure 6. Luminosity, compared with experimental data from Sgr A*, in the case of 
the spiky density profile (left panel) and NFW density profile (right panel). Dashed 
curves represent the luminosity obtained neglecting the effect of advection. The 
computation is performed with = 100 GeV (upper curves) and rriy^ = 1 TeV 
(lower curves). The annihilation cross section is {av)ann = 10"^'' cm^/s. 



SSA decreases the emitted radiation in the frequency range 10^° — 10^^ Hz. This effect 
becomes relevant only when advection is taken into account. 

In order to constrain the shape of the dark matter density profile at the GC or the 
parameters of the dark matter particle, we compare the calculated luminosities with the 
observations. The well-studied source at the GC, Sgr A*, has been observed in a large 
range of frequencies, from the radio to the near-infrared. Data from 10^ Hz to 10^^ Hz 
PH] are displayed in Fig. IHl In the same figure we show the two cases of spiky (left 
panel) and NFW (right panel) density profiles, for = 100 GeV (upper curves) and 

= 1 TeV (lower curves) and {(Jv)ann = 10^^'' cm^/s. In Fig. IHl we also show the 
results for the cases in which only synchrotron losses are taken into account (dashed 
curves) . 

From Fig. IHl it is clear that a spike in the GC induces much stronger emission 
than the observed flux. Therefore, either there is no spike in the dark matter profile 
or neutralinos are not the dark matter. This conclusion qualitatively confirms what 
found in previous studies P, although the magnetic field and radiative transfer here 
are significantly different from what was assumed in [Oj, which results in numerically 
different conclusions (see next section for a more detailed discussion of the comparison 
of our results with those obtained in previous literature). 

The parameter space for neutralinos (m^, {av)ann) has been studied extensively 
jHl EHj. From accelerator and cosmological constraints niy^ G [37 GeV, 7 TeV] and 
(c"'y)ann~3 X 10~^^cm^/s. Siuce the luminosity scales with the neutralino mass and the 
annihilation cross section through the electron density, and n{E,r) oc {av)ann m~^^'^, 
the luminosity we calculated can be lowered by more than two orders of magnitude 
when is increased from 100 GeV to 7 TeV. The luminosity increases linearly for a 
fixed as {av)ann increases. A factor of 30 increase is within the allowed parameters 
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as {(Jv)ann IS changed from 10^^'' cm^/s to 3 x 10"^^ cm^/s. Exploring the parameter 
space can change our results from the fiducial choices (100 GeV, 10~^^ cm^/s) by two 
orders of magnitude in either direction. 

The large enhancement of the neutralino density due to the spiky profile produces 
a synchrotron luminosity that is difficult to reconcile with observations. Moreover, this 
conclusion only gets stronger when advection is included. The SSA effect does not affect 
the discrepancy between a neutralino spike and the observations, since it changes the 
luminosity only at frequencies in the range lO^'' — 10^^ Hz. 

The situation is completely different for the case of the NFW density profile. In 
this case, the synchrotron luminosity is always less than the experimental data as can 
be seen in Fig. NFW is consistent with observations even if {cFv)ann = 3 x 10~^^ 
cm^/s is considered. 

In this paper, we only considered two possibilities for the dark matter density 
profile: a less concentrated hypothesis (NFW) and an extremely concentrated central 
region (spike). As we stressed above, the spike we considered results from the adiabatic 
compression of a NFW dark matter profile in the gravitational field of the central BH. 
Less extreme spikes can be generated if the dark matter is initially distributed in a 
shallower manner compared with the NFW profile. On the other hand, there are cuspier 
original profiles, such as the Moore et al. profile (which scales as r~^/^). A simple Moore 
profile should generate a central luminosity in between the ones illustrated in this paper, 
while a central BH in such a profile will result in more pronounced spike. In some of 
these alternative cases, observations are likely to place more stringent limits on the 
neutralino parameters. 

8. Comparison with previous results 

In this section, we reconsider the case of a spike with magnetic fields as in j^l IT^ . 
We do this to show that these previous calculations have been carried out neglecting 
important pieces of physics, that have now been introduced. We recall that the magnetic 
field assumed in [HI Ej is much smaller than the real equipartition field defined in our 
Eq. Q, although it was quoted in these previous papers as the equipartition field. This 
lower field could in principle be obtained assuming considerably lower accretion rate 
than is usually done. 

The results of jHl 1^ at high frequency are obtained ignoring the effect of advection 
and taking only synchrotron emission as the main channel for energy losses of the 
electrons with a magnetic field as in Eq. (fTUj) . However, as we already discussed in 
previous sections, for such low fields, the effect of advection compared with synchrotron 
energy losses becomes very important, as follows from Eq. ((201): ignoring the 
contribution of advection results in the incorrect prediction for the luminosity of the 
generated radiations. 

Another important point in the comparison with the results of 1^1 is related to 
the role of Synchrotron Self-Compton scattering, which was also neglected in [HI Ej- 
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Figure 7. Comparison between the energy density of the synchrotron emitted photons 
Uph (soHd Hne) and the energy density associated to the magnetic field -B^/Stt (dashed 
line), obtained fixing Bq = 1 fiG and ignoring the effect of advection (as in piT^l. 

While SSC can indeed be neglected when using the correct equipartition magnetic field, 
as discussed in sectional (see Fig. Oj), this is not true in the case of a lower magnetic 
field (cfr. Eq. firUj) ) as the one used in [Hlllll- In this latter case, SSC plays a crucial 
role. This can be easily shown by comparing the energy density of synchrotron photons 
Uphir) with the magnetic energy density Bf^^/Sn. The result is shown in Fig. [3 obtained 
adopting = 10^ GeV and {av)ann = 10"^'' cm^/s. From this figure it is evident that 
the effect of SSC is not negligible at distances larger than ~ 10^/2g away from the GC. 
The effect of SSC is to lower the electron density outside about 10^ Rg from the GC. 

We confirm that the spectrum of synchrotron radiation is affected by strong 
synchrotron self-absorption at the low frequency end, as also pointed out in [T2]. This 
is an important point since in [HI SSA was not included, which led to an overestimate 
of the low frequency synchrotron luminosity. 

In Fig. IHl we neglect advection and compare observations with the emission as 
predicted including SSC, for = 10^ GeV and {av)ann = 10~^^cm^/s. The effect of 
SSC is not enough to reconcile the predicted luminosity with observations. The spike 
hypothesis remains inconsistent with observations. 

9. Conclusions 

In conclusion, we reaffirm that a spike density profile obtained from the adiabatic 
compression of a NFW profile in the gravitational field of the BH at the GC, using dark 
matter parameters typical of neutralinos, is ruled out by the radio and NIR observations 
from SgA* while the NFW density profile remains a suitable hypothesis, not strongly 
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Figure 8. Luminosity obtained ignoring the effect of advcction and taking into 
account the effect of SSC (continuos hne) or ignoring it (dashed hne), compared with 
experimental data from Sgr A* , in the case of the spike density profile with Bq = 1 
liG. 



constrained by available observations. 

We reached this conclusion by a careful consideration of the accretion flow and 
the loss processes in the transport equations for the neutralino generated electrons 
and positrons as well as the radiative transfer. At the Galactic Center, the accretion 
flow onto the central BH sustains strong magnetic fields that induce synchrotron 
emission by electrons and positrons generated in neutralino annihilations during 
advection onto the BH. We found that the magnetic fields associated with the flow 
are significantly stronger than previously adopted. With these equipartition fields, 
we derive the appropriate equilibrium distribution of electrons and positrons and the 
resulting radiation considering adiabatic compression in the accretion flow, inverse 
Compton scattering off synchrotron photons (synchrotron self-Compton scattering), 
and synchrotron self-absorption of the emitted radiation. We calculate the signal for 
a NFW dark matter halo profile and a NFW profile with a dark matter spike due to 
the central BH. We found that the annihilation of neutralino dark matter in the GC 
results in radio signals that overwhelm observations if the spike is present in the dark 
matter density profile. However, observed emissions from the GC are consistent with 
neutralinos following a NFW profile at the Galactic Center. 
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